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Discovery of atomic and molecular mid infra-red emission lines in 
off-nuclear regions of NGC 1275 and NGC 4696 with the Spitzer 
Space Telescope 
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ABSTRACT 

We present Spitzer high-resolution spectra of off-nuclear regions in the central cluster galax- 
ies NGC 1275 and NGC 4696 in the Perseus and Centaurus clusters, respectively. Both objects 
are surrounded by extensive optical emission-line filamentary nebulae, bright outer parts of 
which are the targets of our observations. The 10-37jLim spectra show strong pure rotational 
lines from molecular hydrogen revealing a molecular component to the filaments which has 
an excitation temperature of ^ 300 — 400K. The flux in the O-OS(l) molecular hydrogen line 
correlates well with the strength of the optical lines, having about 3 per cent of the HaH-[NII] 
emission. The 11.3/xm PAH feature is seen in some spectra. Emission is also seen from both 
low and high ionization fine structure lines. Molecular hydrogen cooler than ^400K domi- 
nates the mass of the outer filaments; the nebulae are predominantly molecular 

Key words: galaxies: clusters: general - galaxies: clusters: individual: NGC 1275 - galaxies: 
clusters: individual: NGC 4696 — intergalactic medium - infrared: galaxies 



1 INTRODUCTION 



The massive central galaxy in many clusters i s often surrounded 
by an extensive emi s sion-line nebulosity ( e.g. Cowie et al 19831 : 
I Johnstone et"ail 1 1981 iHeckman et^ll989l : ICrawford et allll999h . 



Such emission-line galaxies are always at the centre of highly 
peaked X-ray emission from an intracluster medium where the ra- 
diative cooling time of the intracluster gas at the centre is less than 
1 Gyr ( lFabianill994l : IPeres et alll 19981 : ISauer et alll2005l) . Initially 
the nebulosity was studied in the optical spectral region but strong 
UV and IR lines have since been seen, together with dust which is 
inferred from the depletion of calcium, Balmer line ratios and dust 
lanes which are prese nt in some cases. Rec e ntly, strong H2 H nes 
have been found (e.g. Ijaffe & Bremed[l997l : lEdge et alll2002l and 
refs therein) as well as emission lines from the CO mole cule (e.g. 
lEdgell200ll ; ISalome & Combes ll2003l : [Salome et alll2006h . 

The ionization, excitation, origin and fate of the filaments has 
been a long-standing puzzle. Most obvious sources of ionization, 
such as an active nucleus, a radio source, shocks, young stars and 
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conduction from the surrounding hot gas have drawbacks in one 
system or another. 



XMM/RGS and Chandra X-ray spectra show that the intra- 
cluster gas in many cluster cores has a range of temperatures de- 
creasing to about a factor of three below that of the outer hot gas, 
with very little X-ray emitting gas cooler than this jPeterson et all 
l200lh . indic ating that some form of heating (e.g. from the central 
radio source. Ichurazov eta l 2002; Fabian et al 2003a) is important 
and that it mostly balances the radiative cooling. Heating is unlikely 
to completely balance cooling over such an extended region, and so 
some residual cooling flow is expected to occur Optically emitting 
warm gas filaments may be being dragged out from the central clus- 
ter galaxy by buoyant bubbles (Hatch et al 2006) or may be a phase 
through which some gas passes in cooling out from the intraclus- 
ter medium. The wider relevance of this issue is in the upper mass 
limit to the total stellar component of massive galaxies, which ma y 
be controlled by whatever heats cooling flows jpabian et all2 002b'). 
If radiative cooling is not balanced in some way (or the cooled 
gas does not form stars) then the stellar component of cD galax- 
ies should be much more massive than is observed, merely by the 
accretion of cooled intracluster gas. By measuring the quantity of 
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young stars and warm/cold gas in these objects we can determine 
how well heating balances cooling. 

The emission-line filaments are markers of feedback in mas- 
sive galaxies and to use them as such we need to understand their 
composition, origin and lifetime. Here we study two of the bright- 
est and most extensive emission-line nebulae around central cluster 
galaxies NGC 1275 (at redshift z = 0.0176, in the Perseus Clus- 
ter) and NGC 4696 (redshift z = 0.009867, in the Centaurus clus- 
ter) using the Spitzer Infrared Spectrograph (IRS). The Ha fila- 
ments around NGC 1275, which extend well beyond the body of the 
galaxy, have warm (2000K) molecular hydrogen associated with 
them. This motivates us to search for much cooler (few x lOOK) 
molecular hydrogen in these regions using the rotational emission 
lines accessible to the IRS. 

The central cD galaxy in the nearby Perseus cluster, 
NGC 1275, is surrounded by spectacular Ha fila ments which 
stretch over 100 kpc. These were first reported by Minkowski 
( Il957l) and imaged bv iLyn ds (1970), McNamara et al ( 1996) and 
recently bv lConselice et ail feoOlh . The Perseus cluster is the bright- 
est cluster of galaxies (in flux uni ts) at X-ray wavelengths with th e 
emission peaking on NGC 1275 jpabian et allll98lLl200dl2003ah . 
The nucleus of NGC 1275 also powers an FRI radio source, 3C 84 
jPedla r et al 1990). Bubbles of relativistic plasma have been blown 
by the jets from the nucleus and have displaced the X-ray emit- 
ting intrac luster gas teohring er et all 1 19931 : IChurazov et a i l200ll : 
iFabian et al 2000, 2002a, 2003 J. It is likely that the Ha filaments 
avoid the bubbles in three dimensional space, although there is no 
obvious anti-correlation as seen in projection on the sky; They ap- 
pear to act like streamlines showing that the flow behind the outer 
bubble, and thus the whole inner medium, is not turbulent and may 
well have bee n dragged out from the central region ( Fabian et al 



2003b, Hatch et al 2006). Molecular hydrogen (Krabbe et al 2000; 
Jaffe et aL2001, ; ,Edge et aL2002) and CO tSalome et al2006, and 



references therein) have been observed and mapped in the inner 
regions for some time. 

NGC 4696 at the centre of th e Centaurus cluster has Ha fil- 
aments discovered by Fabian et all f l982) and a strong dust lane 
llShobbrooklll966t) . ISparks et all l ll989.) and .Crawford et all iimj) 
have mapped the nebulosity which appears to surround the dust 
lane. 

The ionization state of the Ha filaments in central cluster 
galaxies is low, with [OI]A6300, [OII]A3727 and [NII]16584 be- 
ing prominent in optical spe ctra as well as Balmer lines. Molecu- 
lar H 2 is also common (e.g. |jaffe&Breme 3 ll99l I Donahue et al 
2000; Edge et al 2002), even in the outer filaments ( Hatch etal 
2005) . Egami et al ( 2004) have recently found very strong rota- 
tional H2 lines from the brightest cluster galaxy in Zw3146, us- 
ing Spitzer IRS data. Unlike with our objects, the emission is not 
spatially resolved due to the high redshift of the object (z = 0.29). 
The authors remark that this source has the most luminous pure 
rotational H2 lines and the largest mass of molecular hydrogen 
known in a brightest cluster galaxy. A S pitzer IRS spec t rum o f 
the centre of NGC 1275 is presented bv IWeedman et all ( |2005|) . 
Forbidden atomic lines and at least one molecular line are seen 
sitting on a l arge continu um associated with the active nucleus. 
iKaneda et al | |2005| , l2007h have taken low resolution Spitzer IRS 
spectra of NGC 4696 and find PAH features and [Nell] line emis- 
sion. 

What is ionizing and exciting the fila mentary gas remains 
uncertain. It is not sim ply an active nucleus ( jjohnstone & Fabiat] 
ll988l ; ISabraetaill2000l) . nor do a range of alternative mechansims 
provide a comprehensive explanation (e.g. ICrawford & FabiMil 



I1992I ; jPonahue et aill2000l ; [ Sabra et a 3 120001 ; IWilman et alll2002h . 

The UV emissi on from massive young stars is likely to play an 
impo rtant role ( jjohnstone et alll 19871 ; lAllenlll99^ ; ICrawford eTail 
1 19991) , although this may not be the case for the outer filaments. 
The lifetime of the filaments is not known, although the immense 
size of the system around NGC 1275 (80 kpc) and the low velocity 
spread seen in the filaments (~ 300 km s^' or less) indicates that 
they may last over 10^ yr. 

In this paper we present mid infra-red Spitzer spectra of 
off-nuclear extended emission-line regions in NGC 1275 and 
NGC 4696. They probe a new region in temperature space (~ 
300 — 400K) which lies between the very cold regions observed 
in transitions of the CO molecule and the warmer regions seen 
though ro-vibrational transitions of the hydrogen molecule. These 
emission-line filament s seem to be a key marker of feedb ack in 
cooling core clusters dHu et all [1983; ICrawford et all 1 19991) even 
though they do not contain a significant fraction of the total gas 
mass. It is by studying regions located away from the nucleus of 
the galaxies that we expect to be able to make progress in the un- 
derstanding of the heating and ionization mechanisms of these fil- 
aments as they offer a simpler environment than is found close to 
the nucleii of the galaxies. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Mid-infrared Spitzer Data 

Spitzer IRS observations of off-nuclear regions in NGC 1275 and 
NGC 4696, were made during two observing campaigns in 2005 
July and 2006 March (IRSX005200 and IRSX006300). For each 
object both the Short Wavelength High Resolution (SH) and Long 
Wavelength High Resolution (LH) spectrographs were used, giving 
an effective resolving power of R~ 600 and a spectral coverage 
from ~ 10/xm to ~ 37;tim. The entrance aperture size of the two 
spectrographs is rather different, being 4.7 x 1 1 .3 arcsec for the SH 
and 1 1 . 1 X 22.3 arcsec for the LH spectrographs respectively. 

Table [T] lists details of the Spitzer observations. Three off- 
nuclear regions were observed in NGC 1275, chosen to be coin- 
cident with particular regions identified in the Ha-(-[NII] map of 
IConseliceet"ail (l200l'). as well as a blank sky region. In NGC 4696 
we covered the brightest part of the Ha nebula to the South West 
of the nucleus. 

Figures [T] and [2] show the outline of the Spitzer spectrograph 
apertures overlaid on Ha-l-[NII] images of the central cluster galax- 
ies. There are two slightly different pointings (exposure ids) for 
each target position, these being the standard nod positions in the 
Staring Mode Astronomical Observation Template. These positions 
locate the target at one third and two thirds of the way along the 
spectrograph aperture. The SH spectrograph aperture positions are 
shown by the smaller rectangles whereas the LH spectrograph aper- 
tures are the larger rectangles which are approximately orthogonal 
to the SH spectrograph apertures. The specific pointing positions 
for each exposure id are given in table[T] 

The observations were processsed by the Spitzer Science 
Center through pipeline version S13.2.0. We use the 'bcd.fits', 
'func. fits' and 'bmask.fits' files as the starting point for our reduc- 
tion. 

Individual Data Collection Event (DCE) files within each ex- 
posure id were first averaged together and an array of external un- 
certainties corresponding to the standard error on the mean was 
calculated. A new bmask file was also generated, this being the log- 
ical 'or' of all the bmasks for the input DCEs. The NGC 1275 East 
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Table 1. Log of Spitzer Observations. 



Target Redsliift Obs date AOR key Mode Exposure id RA Dec Exposure time 

(2000) (seconds) 



JNLtL. Iz /j nasi 


0.0176 


2UUo Mar 14 


.-1/1 ^1 ^,1 nA 
rl4j JO /04 


SH 


0002 


03: 


19:49.9 


1 A 1 

+41 


:30:48 


481.69 












0003 


03: 


19:50.2 


+41 


:30:46 


481.69 










LH 


0004 


03: 


19:49.9 


1 A 1 

+41 


• In- A A 


241.83 












0005 


03: 


19:50.2 


+41 


:30:50 


241.83 


JNCjC 12/5 Position 2 


0.0176 


2006 Mar 14 


rl45io44o 


SH 


0002 


03: 


19:45.7 


+41 


:29:49 


1445.07 












0003 


03: 


19:46.0 


+41 


:29:47 


1445.07 










LH 


0004 


03: 


19:45.7 


+41 


:29:45 


725.48 












0005 


03: 


19:46.0 


+41 


:29:51 


725.48 


NGC 1275 Position 11 


0.0176 


2006 Mar 14 


rl4536192 


SH 


0002 


03: 


19:45.0 


+41 


:31:34 


1445.07 












0003 


03: 


19:45.3 


+41 


:31:32 


1445.07 










LH 


0004 


03: 


19:45.0 


+41 


:31:30 


725.48 












0005 


03: 


19:45.3 


+41 


:31:36 


725.48 


NGC 1275 Blank 


0.0176 


2006 Mar 14 


rl6754176 


SH 


0000 


03: 


19:53.6 


+41 


:28:35 


963.38 












0001 


03: 


19:53.9 


+41 


:28:33 


963.38 










LH 


0002 


03: 


19:53.6 


+41 


:28:31 


483.65 












0003 


03: 


19:53.9 


+41 


:28:37 


483.65 


NGC 4696 


0.009867 


2005 Jul 9 


rl4537216 


SH 


0002 


12:48:48.7 


-41: 


18:42 


2408.44 












0003 


12:48:48.9 


-41: 


18:45 


2408.44 










LH 


0004 


12:48:48.5 


-41: 


18:45 


725.48 












0005 


12:48:49.1 


-41: 


18:42 


725.48 




Figure 2. Spitzer spectrograph a pertures overlaid on a n Ho;+[Nn] image of 
Figure 1. Spitzer spectrograph apertures overlaid on an Ha+[NII] image of (he central region of NGC 4696 jcrawford et all2"ooi) . North is up and East 

the central region of NGC 1275 f Conselice et a<|200l | ) . Small boxes show ^^^^ j^ft g^^u ^oxes show the positions of the SH spectrograph aperture 

the positions of the SH spectrograph aperture while the larger boxes show ^yj^ ^^^^^ boxes show the positions of the LH spectrograph aperture, 

the positions of the LH spectrograph aperture. 



pointing has only one DCE in each spectrograph so in these cases 
the pipline based uncertainties, which are computed from the ramp 
fitting of the multiple detector readouts, are used for the uncertainty 
arrays. 

The NGC 1275 pointings were then background subtracted us- 
ing the NGC 1275 blank field pointing, propagating the uncertainty 
arrays in quadrature and again creating a new bmask file that is the 
logical 'or' of the input masks. 

In the case of the NGC 4696 pointing we do not have a con- 
current sky-background observation. For the LH data, we used an 
average of background observations (scaled to the expected back- 
ground level at the position of the NGC 4696 pointing using val- 
ues from the Spitzer SPOT software) taken for standard star cali- 



brations in the same observing campaign to effect sky subtraction. 
However, for the SH data, the analagous background observations 
were too short compared with the NGC 4696 data and added con- 
siderable noise when subtracted. We therefore leave these as not 
sky subtracted. One consequence of not having a concurrent sky 
background observation is that the rogue pixels in the detectors are 
not well corrected (or corrected at all in the case of the SH data). 
We therefore used the IRSclean_mask software from the Spitzer 
contributed software sit^ to interpolate over the rogue pixels for 
these observations. 

The data from the SH and LH spectrographs are not rectilinear 



' |http : //ssc ■ spitzer . caltech . edu/archanaly/ contributed/browse ■ htr 
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Table 2. Ratio of H g+pSHI] fluxes i n the Spitzer apertures measured from 
the or iginal maps of lConselice et alH200ll) (NGC 1275) and lCrawford et all 
(2005) (NGC 4696) to fliose smooflied with a gaussian kernel of FWHM of 
3.10 arcsec (for the SH spectra) and a FWHM of 6.17 arcsec (for the LH 
spectra). 



Region 


Exposure id 


Spectrograph 


Flux ratio 


JNtjC 111 J xiast 


0002 


SH 


1.00 




UUUj 




i.Uo 




0004 


LH 


uo 




0005 


LH 


Ul 


NGC 1275 Position 2 


0002 


SH 


1.07 




0003 


SH 


1.13 




0004 


LH 


1.04 




0005 


LH 


1.08 


NGC 1275 Position 11 


0002 


SH 


1.09 




0003 


SH 


1.11 




0004 


LH 


1.07 




0005 


LH 


1.01 


Centaurus 


0002 


SH 


0.99 




0003 


SH 


1.00 




0004 


LH 


1.16 




0005 


LH 


1.06 



in either the spatial or spectral direction and therefore require soft- 
ware which has a knowledge of the detailed distortions in order to 
be able to extract one-dimesional spectra. We used the Spitzer Spice 
software (version 1.3betal and version 1.4) for this, extracting all 
the spectral orders over the full aperture in both spectrographs. The 
final step in the Spice software extraction is a 'tuning' that applies a 
flux calibration to the data. There are two possibilities at this stage 
because the size of the point-spread function of the Spitzer tele- 
scope is of the same order as the entrance aperture size. One option 
applies an extended source calibration in which it is assumed that 
the aperture is illuminated by a uniform surface brightness source; 
as much flux is scattered into the aperture as is scattered out of it. 
The other option is to apply the point-source calibration which cor- 
rects for the fraction of flux from a well-centered point-source that 
is lost from the aperture, as a function of wavelength. 

To determine whether the point-source or extended source cal- 
ibration is more appropriate for our data we have taken the contin- 
uum subtracted Ho!-l-[NII] maps and used them as proxies for the 
spatial distribution of the infra-red line emission since we do not 
know its spatial distribution. Assuming that the Ho!-l-[NII] emis- 
sion is distributed spatially like the infra-red emission allows us 
to make a quantitative (but approximate) assessment of whether 
the emission, which appears very clearly extended on the resolu- 
tion of the optical data, is actually extended at the resolution of the 
Spitzer observations. We have measured the ratio of the flux in the 
Spitzer apertures in the raw image and in the image smoothed using 
a gaussian kernel of full- width half-maximum of 3. 1 or 6. 17 arcsec. 
These widths are an approximation to the size of the Spitzer point 
spread function at 15/im and 30/xm (approximately the mid-points 
of the SH and LH spectra respectively) as measured by fitting a 
one-dimensional gaussian model to the azimuthal average of the 
central peak of a point source generated by the Spitzer Stinytim 
point-spread function modelling software. Table|2]lists these ratios 
for each of the exposure ids. 

Examination of Table [2] shows that there is up to 16 per cent 
more Ha-l-[NII] flux measured in the Spitzer apertures in the un- 
smoothed images compared with the smoothed images. The in- 



verse of the Spitzer Slitloss Correction Functior0 shows that for 
a point source observed by Spitzer the factor would be 1.63 and 
1.52 at 15/im and 30/lm in the SH and LH spectrographs respec- 
tively, whereas a uniform brightness extended source would give 
factors of 1.0 for both wavelengths. This indicates that the data ap- 
proximate much better to an extended source than a point source. 
We therefore proceed using the extended source tuning option in 
Spice. 

The final part of the reduction involves merging the individ- 
ual spectral orders into one spectrum and converting the flux units 
to per unit wavelength interval. The sensitivity of the spectrograph 
drops significantly at the end of each order leading to an increase 
in the noise, but because there is an overlap in the wavelength cov- 
erage of adjacent orders (and the wavelength binning is the same at 
corresponding wavelengths) we are able to use the following pre- 
scription to clean the final spectrum: In the order overlap regions 
we choose data from the lower order number spectrum in prefer- 
ence to that of the higher order number, except that we discard the 
6 shortest wavelength bins in the highest order spectrum and the 
4 shortest wavelength bins of all other orders. The 6 longest wave- 
length bins in the spectrum are also discarded. The spectral binning 
used by the Spice extraction routine varies as a function of wave- 
length to preseve a constant velocity resolution of approximately 
240 km s^' per bin. We show examples of the final spectra in Figs. 
(3) |4l |5] and |6l with the expected positions of some emission lines 
marked. The dashed vertical lines show the boundary of adjacent 
spectral orders. 

The emission lines in the spectra were fitted individually using 
a gaussian model for the line plus a straight line to account for the 
local continuum. The fitting was done using the QDP prograrr0 to 
minimize the chi-square statistic between the model and the data, 
subject to the propagated uncertainty in the flux at each wavelength 
bin. Line parameters and uncertainties are given in Tables [3] ID [5] 
and|6] The uncertainties were generally calculated (for all pointings 
except the NGC 127 5 East region) from the l^x^ = 1.0 criterion 
jLampton et all 197^ . This corresponds to a 68 per cent confidence 
region or a IcT (gaussian equivalent) confidence region for one in- 
teresting parameter. In the cases where lines are not detected we set 
a 30^ upper limit on the line flux by fixing the position and width 
of the line and increasing the flux until the value of chi-square in- 
creases by 9 from the value obtained with no line present. Where 
no uncertainty on a parameter is listed that parameter was fixed. 
Line widths for detected lines were left as free parameters in the 
fitting. However, except for the PAH I1.3/xm line the widths are 
not resolved beyond the instrumental resolution and are therefore 
not listed. 

The values of the minimum chi-square that were obtained 
when fitting the NGC 4696 short wavelength spectra were signif- 
icantly above the number of degrees of freedom, whereas the val- 
ues of the minimum chi-square that were obtained when fitting the 
NGC 1275 East spectra were significantly below the number of de- 
grees of freedom. In both cases this indicates that the uncertainties 
on each bin are wrong. In the case of the NGC 4696 short wave- 
length data there is additional power in the spectra not accounted 
for in the uncertainty arrays which may be due to the presence of 
fringes (http://ssc.spitzer.caltech.edu/postbcd/irsfringe.html ) that 
have not been completely removed in the pipleline sofware. We 
tried to correct for these using the irsfringe tool but were unsuc- 



^ http : // ssc . spit zer . calt ech . edu/ irs/calib/ ex tended-sources /index 
^ http : //wwwastro .msfc.nasa. gov/qdp 
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Figure 3. Upper: Spitzer SH (exposure id 0002) and Lower: Spitzer LH 
(exposure id 0004) spectra for the NGC 1275 East pointing. Daslied vertical 
lines show the boundary between adjacent spectral orders. 
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25 30 
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Figure 4. Upper: Spitzer SH (exposure id 0002) and Lower: Spitzer LH 
(exposure id 0004) spectra for the NGC 1275 Position 2. Dashed vertical 
lines show the boundary between adjacent spectral orders. 



cessful. In the case of the NGC 1275 East spectra there is only one 
DCE and the uncertainties are propagated from the ramp fitting. 
This leads to uncertainties which are too big. 

Since using the A;^^ = 1.0 criterion under either of these con- 
ditions leads to uncertainties on the line parameters which are ei- 
ther too la rge or too small we h ave adopted the 'Ratio of Variances' 
technique jLampton et a lll976h to rescale the critical value of delta 
chi-square used to calculate the line parameter uncertainties. We 
note in passing that the spectral extraction method used in the Spice 
software does introduce correlations between adjacent spectral bins 
but the Spitzer Science Centre advises that the degreee of correla- 
tion is small because the dispersion direction in the spectral images 
is almost parallel with the pixel direction. 



2.2 Supporting Near-infrared Data 

The Spitzer apertures at positions 2 and 11 in the NGC 1275 neb- 
ula cover similar sky positions to the A^— band United Kingdom 
Infrared Telescope (UKIRT) CGS4 long - slit observations of the 
'horseshoe knot' and the 'SWl' region of iHatch et"^ j2005l) . Not- 
ing that the entrance apertures for the Spitzer spectrographs are a 
different shape and size to the region extracted from the CGS4 slit 
we will later use the Po! and H2 ro -vibrational line fluxes of these 
two regions from lHatch et al (l2005h together with the Spitzer mid- 
infrared H2 rotational line fluxes of positions 2 and 1 1 to investi- 
gate the molecular hydrogen in these outer regions of the nebula. 



Full details of the observations and d ata reduction of th e A'-band 
data from these regions are provided in lHatch et all ( l2005h . 

The Spitzer aperture covering the East reg ion does not match 
well to the Eastern slit region of iHatch et all 12005) . Instead we 
have used UKIRT UIST //A'-band longslit spectra of the bright 
radial filament that extends from a region 2.65 arcsec South of 
the galaxy nucleus eastward to a distance of 26.8 arcsec (9.5 kpc). 
The observations were obtained on 2005 January 31, and the 
sky conditions were photometric. The long-slit was used in a 
7 pixel-wide (0.84 arcsec) configuration with a spatial scale of 
0.1202 arc sec pixel ^ using the ///f grating which gives a spectral 
coverage of 1.4 — 2.5/xm at a resolution of 700— 1260 km s^' . The 
observations were taken in NDSTARE mode with an object-sky- 
sky-object nodding pattern. Details of the observation are provided 
in TablelD Sky emission features were mostly removed by the nod- 
ding pattern, while flux calibration was achieved using a set of al- 
most featureless F and G type stars of known magnitude. In each of 
the standards a Bry stellar feature at 2. 166/xm was removed by lin- 
ear interpolation. The data were reduced with ORAC-DR available 
through the UKIRT Website, and Starlink packages. 

A spectrum was extracted from the UIST long-slit data by 
summing the flux across the entire region where the more northerly 
Spitzer aperture overlaps with the UIST aperture. This results in 
a total long-slit aperture of 0.84x6.6 arcsec^. Fig. [7] shows the 



' |http : / /www ■ jach . hawaii ■ edu/UKIRT/ instruments/uist/ spectroscopy/ 
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Table 3. Properties of emission lines in the NGC 1275 East region Spitzer spectra. Note that the line widths are the gaussian sigma values. 



Exposure id 



Line 



Wavelength 
(microns) 



Width (&) 
(microns) 



(xlO erg cm 



Flux 

= s-') 



0002 



0003 



0004 



0005 



PAH 11.3 
H2 0-0S(2) 12.28 
[Nell] 12.81 
[Nelll] 15.56 
H2 O-OS(l) 17.04 
[Sill] 18.71 
PAH 11.3 
H2 0-0S(2) 12.28 
[Nell] 12.81 
[Nelll] 15.56 
H2 O-OS(l) 17.04 
[Sill] 18.71 
H2 O-OS(O) 28.22 
[Sill] 33.48 
[Sill] 34.82 
H2 O-OS(O) 28.22 
[SlII] 33.48 
[Sill] 34.82 



11.51±0.01 
12.4899±0.0003 
13.0342±0.0002 
15.8241±0.0007 
17.3280±0.0009 
19.040±0.003 
11.56±0.03 
12.4906±0.0003 
13.0350±0.0003 
15.8236±0.0010 
17.3270±0.0010 
19.()40±0.003 
28.72 
34.07 
35.421±0.001 
28.72 
34.07 
35.415±0.001 



O.lliO.Ol 



0.14±0.03 



1.43±0.16 
0.85±0.03 
1.78±0.03 
0.54±0.03 
1.79±0.()9 
0.39±0.07 
1.15±0.30 
0.87±0.03 
1.51 ±0.04 
0.44±0.03 
1.98±0.12 
0.16±0.04 
<0.73 
<1.68 
4.46±0.18 
<0.75 
<0.75 
3.59±0.17 



Table 4. Properties of emission lines in the NGC 1275 Position 2 region Spitzer spectra. Note that the line widths are the gaussian sigma values. 



Exposure id 


Line 


Wavelength 


Width (a) 


Flux 






(microns) 


(microns) 


(x 10"'* erg cm"- s"' ) 


0002 


PAH 11.3 


11.51 


0.11 


<0.11 




H2 0-0S(2) 12.28 


12.4920±0.0008 




0.24±0.02 




[Nell] 12.81 


13.0358±0.0006 




0.30±0.02 




[Nelll] 15.56 


15.8241 




<0.19 




H2 O-OS(l) 17.04 


17.33 17±0.0004 




0.47±0.02 




[SlII] 18.71 


19.040 




<0.10 


0003 


PAH 11.3 


11.51 


0.11 


<0.11 




H2 0-0S(2) 12.28 


12.4925±0.0010 




0.27±0.03 




[Nell] 12.81 


13.0382±0.0011 




0.21±0.03 




[Nelll] 15.56 


15.8241 




<0.11 




H2 O-OS(l) 17.04 


17.3334±0.0006 




0.42±0.02 




[Sill] 18.71 


19.04 




<0.07 


0004 


H2 O-OS(O) 28.22 


28.72 




<0.15 




[SIII] 33.48 


34.07 




<0.27 




[Sill] 34.82 


35.43 




<0.63 


0005 


H2 O-OS(O) 28.22 


28.72 




<0.18 




[SIII] 33.48 


34.07 




<0.36 




[Sill] 34.82 


35.43 




<0.69 



Table 8. Details of the UIST long-slit observation of the NGC 1275 East 
region. The orientation of the nodding direction is described in the UIST 
manual' . Position angle is measured East from North. 



Target position 


Slit angle 


Nod offset 


Exposure 


Seeing 


(J2000) 


(degrees) 


(arcsec) 


time (mins) 


(arcsec) 


03 19 48.16 


+70 


+38, -23 


144 


0.4 


41 30 40.1 











NGC 1275 East near-infrared spectrum over the wavelength range 
1.74 — 2.48/im. The wavelength range 1.8— 1.88/xm is a region of 
poor atmospheric transmission and has been masked out of the 
spectrum for clarity of the emission-line features. The molecular 
hydrogen emission lines and Pa are clearly visible. Table |9] lists 



the surface brightnesses of all the emission-lines. These are lower 
limits to the intrinsic surface brightnesses because of the unknown 
covering factor of the excited molecular hydrogen gas. 



3 ANALYSIS 

The lines detected in the Spitzer mid-infrared spectra consist of 
emission from molecular hydrogen, from various atomic species of 
Neon, Sulphur and Silicon and from the 11.3 micron PAH featiu-e. 

3.1 Diagnostic line ratios 

In order to classify the emission-line regions seen by Spitzer, we 
have plotted the lines of N eon/Sulphur/Silic on in Fig [8] on one of 
the diagnostic diagrams of lOale etall ( |2006|) . The NGC 1275 East 
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Table 5. Properties of emission lines in the NGC 1275 Position 1 1 region Spitzer spectra. Note that the line widths are the gaussian sigma values. 



Exposure id 


Line 


Wavelength 


Width (a) 


Flux 






(microns) 


(microns) 


(x 10"^"^ ergcm"^ s~^ ) 


UUUz 


PA tr 1 1 T 
Jr7\rl I L.J 


i 1 . J i 


n 1 1 

U. 1 i 


< u.z / 




rl2 U-Uo(^ZJ iZ.Zo 


iZ.^VJJZtU.UUi J 




U. i JitU.UZ 




[INellJ iz.fSi 


1 Q n/im_i_n nnnc 




u.zyihu.uj 




[INelllJ ij.jo 






<U.U / 




rl2 U-Uo(l ) i /.U4 


i /. JJ4z±U.UUUo 




U.zjihU.Uz 






1 o n/1 n 
iy.U4U 




<U.iD 


0003 


PAH 1 1.3 


1 1.5 1 


0. 1 1 


< u.iy 




rl2 U-Uo(z) IZ..Z.O 


iZ.4yjo±U.UUiZ 




A 1 1_LA AT 

U.i J±U.Uz 




[Nell] 12.81 


13.0411±0.0011 




0.15±0.02 




[Nelll] 15.56 


15.8241 




<0.08 




H2 O-OS(l) 17.04 


17.3371±0.012 




0.21±0.02 




[SIII] 18.71 


19.04 




<0.11 


0004 


H2 O-OS(O) 28.22 


28.72 




<0.07 




[SIII] 33.48 


34.07 




<0.38 




[Sill] 34.82 


35.43 




<0.49 


0005 


H2 O-OS(O) 28.22 


28.72 




<0.14 




[SIII] 33.48 


34.07 




<0.24 




[Sill] 34.82 


35.43 




<0.20 



Table 6. Propeilies of emission lines in the NGC 4696 Spitzer spectra. Note that the line widths are the gaussian sigma values. 



Exposure id 


Line 


Wavelength 


Width (fj) 


Flux 






(microns) 


(microns) 


(x 10^'** ergcm"" s^' ) 


0002 


PAH 11.3 


11.32±0.017 


0.15±0.02 


1.44±0.2 




H2 0-0S(2) 12.28 


12.408±0.003 




0.17±0.03 




[Nell] 12.81 


12.9469±0.0003 




1.13±0.03 




[Nelll] 15.56 


15.7184±0.0005 




0.66±0.02 




H2 O-OS(l) 17.04 


17.2122±0.0011 




0.45±0.03 




[SIII] 18.71 


18.91 10±0.004 




0.14±0.03 


0003 


PAH 11.3 


11.33±0.016 


0.11±0.019 


1.75±0.21 




H2 0-0S(2) 12.28 


12.411±0.004 




0.18±0.04 




[Nell] 12.81 


12.9488±0.0003 




0.98±0.03 




[Nelll] 15.56 


15.7207±0.0004 




0.58±0.02 




H2 O-OS(l) 17.04 


17.2142±0.0006 




0.42±0.02 




[SIII] 18.76 


18.9133±0.005 




0.13±0.04 


0004 


H2 O-OS(O) 28.22 


28.497 




<0.12 




[SIII] 33.48 


33.811 




<0.94 




[Sill] 34.82 


35.1808±0.0030 




1.81±0.15 


0005 


H2 O-OS(O) 28.22 


28.497 




<0.10 




[SIII] 33.48 


33.811 




<1.60 




[Sill] 34.82 


35.1740±0.0021 




3.05±0.18 



and NGC 4696 points both lie in the Region I+II area of the di- 
agram occupied by Seyfert and Liner nuclei. This finding is con- 
sistent with other analyses of optical emiss ion lines in these ob- 
jects which als o show Liner-like spectra (eg I Johnstone et alll 19871 : 
ICrawfordetai|[l999i) . 

3.2 [Nelll] lines 

The SH spectra of the NGC 1275 East region and the NGC 4696 
pointing both show strong detections of the [Nelll] Al5.56/xm 
line. In the NGC 1275 East region this is a surprising discovery 
since most of the very extended emission-line gas do es not show 
emission from the optical lines of [OIII]1A4959,5007 faatch et all 
l2006h . The ionization potential of Ne+ (41.07 eV) is larger than the 
ionization potential of 0+ (35.1 leV) so if lines of [Nelll] are seen 
we would expect to also see lines of [OIII]. 



There are two possible explanations for the apparent lack of 
the [OIII] lines when the [Nelll] lines are strong, assuming that 
they are both formed in a photoionized nebula. The first possibility 
is that the [OIII] and [Nelll] lines are both emitted together from 
a spatial region which is missed in our optical long-slit spectra but 
sampled in the Spitzer spectra due to the much larger SH aperture 
(4.7 X 11 .3 arcsec) and the poorer point-spread function of Spitzer. 
Alternatively, since these lines are collisionally excited and the ex- 
citation energy of the upper level in the [OIII] lines (~29,000K) 
is much greater than that for the upper level in the [Nelll] line 
(~940K), if the electron temperature in the emitting gas is low 
enough, the infra-red [Nelll] line may be produced while the op- 
tical [OIII] lines are not excited. 

To determine whether the first possibility is occurring we note 
that the SH aperture for our NGC 1275 E ast observations is lo- 
cated very close to an HII region identified bv lShields & Filippenkd 
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Table?. Ha+[NII] emission-line fluxes in Spitzer apertures measured from maps of lConselice et alH200lh (NGC 1275) an d ICrawford et"3 )2005h (NGC 4696). 
To account for the Spitzer point spread function, the maps were smoothed with a gaussian kernel of FWHM of 3.10 arcsec for the short wavelength aperture 
measurements and by gaussian with a FWHM of 6. 17 arcsec for the long wavelength aperture measurements. No measurement is available for the NGC 4696 
LH apertures as these cover a region of the image in which the off-line image is saturated. The NGC 1275 values have been corrected for a reddening 
conesponding to Ar = 0.4. 



Region 


Exposure id 


Spectrograph 


Flux 

(x 10^'^ ergcm"^ s^' ) 


NLiL iz / J cast 


0002 


SH 


5.8 




UUUj 








0004 


LH 


16 




0005 


LH 


15 


NGC 1275 Position 2 


0002 


SH 


1.1 




0003 


SH 


1.0 




0004 


LH 


2.7 




0005 


LH 


2.8 


NGC 1275 Position 11 


0002 


SH 


0.83 




0003 


SH 


0.68 




0004 


LH 


1.7 




0005 


LH 


1.4 


NGC 4696 


0002 


SH 


1.7 




0003 


SH 


1.5 




0004 


LH 






0005 


LH 






12 14 16 18 

Observed wavelength (Microns) 




25 30 
Observed wavelength (Microns) 

Figure 5. Upper: Spitzer SH (exposure id 0002) and Lower: Spitzer LH 
(exposure id 0004) spectra for the NGC 1275 position 11. Dashed vertical 
lines show the boundary between adjacent spectral orders. 



PAH H, 0-0S(2) 



INellll H,0-0S(1) ISUIl 



12 14 16 

Observed wavelength (Microns) 





20 25 30 35 

Observed wavelength (Microns) 

Figure 6. Upper: Spitzer SH (exposure id 0002) and Lower: Spitzer LH 
(exposure id 0004) spectra for the NGC 4696 pointing. Dashed vertical lines 
show the boundary between adjacent spectral orders. Note that the SH data 
are not background subtracted. 
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Figure 7. Near-infrared spectrum of the NGC 1275 East region. The re- 
gion between 1.8— 1.88^m has poor atmospheric transmission and has been 
masked out of the spectrum for clarity. Errors are from Poisson statistics. 

Table 9. Near-infrared hne surface brightnesses for the NGC 1275 East 
region together with IfJ uncertainties. The value for the 2 — 1 S(l) line is a 
3c7 upper limit. 



Line 




Surface Brightness 

X 10"'* ergcm"^ s"' arcsec"^ 


Pa 




1.73±0.08 


H2I - 


0S(0) 


0.24±0.05 


H2I- 


0S(1) 


1 26+°"** 


H2I- 


0S(2) 


0.55±0.05 


H2I- 


0S(3) 


1.43±0.05 


H2I- 


0S(4) 


0.23±0.08 


H2I- 


0S(7) 




H2I - 


0Q(1) 


1.17±0.11 


H2I- 


0Q(3) 


1.09±0.11 


H22- 


1S(1) 


<0.12 





10.0 

























1.0 




f 


z 


r NGC 4696/ 


1 _ 
1 : 














/ -< • / 




iri 












/ 




z 




NGC 1275 East/ 






0.1 










1/ II ^ II 


IV 



0.1 1.0 10.0 

[SIII]33.48/[Sill]34.82 



Figure 8. Neon/Sulphur/Silicon diagnostic diagram of bale et all i20O6h . 
Regions I-l-II demark the locus of Seyfert and Liner nuclei. Region III de- 
marks the locus of HIT nuclei while Region IV is the locus of Extranuclear 
and HII regions. 



That HII region has much stronger [OIII] emission than the 
surrounding nebula and therefore may contribute significant [Nelll] 
emission. 

In order to assess whether the Spitzer SH apertures contain 
a significant amount of light from the HII region we used the 
Spitzer Stinytim point-spread function modelling software to gen- 
erate an image of the point-spread function for the SH instrument 
at the observed wavelength of [Nelll] (15.8/im). We then applied 
a world coordinate system to this image such that the centre of the 
point spread function lies at the coordinates of the HII region and 
the pixel scale is correct for the point-spread function image. The 
Spitzer aperture region files corresponding to the two exposure ids 
were then applied to this image and the counts in each region were 
measured. 

This analysis predicts that the region corresponding to expo- 
sure id 0002 (the more westerly one, nearer the HII region) should 
have 2.86 times the flux present in the region corresponding to ex- 
posure id 0003 if the flux comes entirely from the point source HII 
region. The ratio of fluxes measured in the [Nelll] lines in these two 
exposure ids is just 1.23. If we further assume that each exposure id 
has a contribution from the extended (non-HII region) nebula pro- 
portional to the Ha-i-[NII] emission listed in Table [T] then we can 
solve for the fractional contribution of the HII region to the [Nelll] 
flux in each exposure id. Under these assumptions, we find that 3 
per cent of the [Nelll] flux (1.5 x 10"'^ erg cm"^ s"' ) in exposure 
id 0002 comes from the HII region while only 1 per cent of the 
[Nelll] flux (5.4 X 10"'^ erg cm"^ s"' ) in exposure id 0003 comes 
from the HII region. The [Nelll] flux from the extended nebula is 
5.25 X 10"'^ ergcm"^ s"' and 4.35 x 10"'^ ergcm"^ s"' in ex- 
posure ids 0002 and 0003 respectively. This suggests that there is a 
strong extended [Nelll] emission component in the NGC 1275 East 
region. 

We have set a 3cT upper limit on the flux in the [OIIIJA5007 
line by extracting a spectrum from the Gemini data, presented by 
iHatch et all feoOO ). close to the Perseus East region (and avoiding 
the HII region). In order to avoid having to correct for different 
spatial region sizes between the optical and infra-red data we pro- 
ceed by taking ratios of the [OIII] and [Nelll] lines to the Ho!-l-[NII] 
lines: 

F(5007) F(5007) F(15.56) 

F(15.56) ^ F(Ha + [Nil]) ' F(Ha + [Nil]) 

where the ratio F(5007)/F(Ha + [Nil]) is the ratio of the 
flux in the [OIII] A 5007 line to the flux in the sum of the 
Ha-i-[NII]A16548,6584 lines both measured from the Gemini 
spectrum. The ratio F(15.56) /F(Ha + [Nil]) is measured from the 
[Nelll] flux in the Spitzer spectrum and the values given in Table 
[3 We find a 3(7 upper limit on the ratio F(5007)/F(15. 56) < 1.2. 

Gas that is photoionized by a stellar or power law 
continuum gene rally has a temperature around 10,000 K 
tos terbrock & Fe rland 2006). The lack of optical [OIII] emission 
at a position where infrared [Nelll] lines are seen suggests that 
the electron temperature is low enough to prevent opt ical emission 
from being produced. Models using the code Cloudy dFerland et all 
Il998l) show that if the abundance of 0+^/Ne+^ has the solar ra- 
tio, then the ratio F(5007)/F(15.56) ~ io4Tj:^/^e(-29000/T,) 
observed limit of 1.2 in this ratio therefore corresponds to Tg < 
6300K. 

Low temperatures in a photoionized gas are produced when 
the abundances are at or above solar and the gas density is low 
enough that most of the cooling luminosity is radiated by th e 
fine structure lines dFerland et allll984IShields & Kennicut^ll995h . 
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These calculations show that metallicities in the range of solar to 
twice solar can produce a drop in the [OIII] line strength by up 
to 2 orders of magnitude depending on the details of the physical 
system. 

The results from the discussion above rely on the observed op- 
tical line ratios (which have been corrected for reddening within our 
Galaxy) not being affected by any obscuration within the emission 
line regions themselves. The Balmer decrement, calculated from 
the ratio of Ha/H/3 in the Gemini spectrum, is very high (6.7) 
even after correcting for Galactic reddening. If this high value of 
the Balmer decrement is attributed to extinction in a screen in front 
of the emission-line region, then the implied extinction is Av=2.67. 
After correcting the optical lines for this amount of reddening we 
obtain a hmit on the ratio F(5007)/F(15.56) < 19 which would 
not set an interesting limit on the electron temperature in the neb- 
ula. We note that high values of the Balmer decrement may be pro- 
duced by mechanisms other than reddening, for example collisional 
excitation (Table 5, Parker 

In NGC 4696, IJohnstone & FabiarJ h988l) report 
[OIII]A5007/Hj3 ~ 0.2 in a 13.2 x 1.5 arcsec region at position an- 
gle 85° centred 5 arcsec south of, but including, the nucleus, while 
iLewis et all ( l2003h set an upper limit of [OIII] A5007/H/3 < 0.4 in a 
region 1.8 x 1.0 arcsec centred on the nucleus. This indicates that 
the [OIII] line is quite weak in this object also. We do not attempt 
to calculate the F(5007)/F(15.56) ratio due to the much larger 
uncertainties because of the proximity of the Spitzer aperture to the 
galactic nucleus (Fig.|2](. The inner region of the Centaurus cluster 
is known to have a hig h metallicity from X-ray observations (eg 
ISanders & Fabiaij|2006l) . although the exact value derived for the 
very central regions is dependent on the precise model fitted. If 
future observations were to be able to constrain the metallicity of 
the emission-line gas, that information might prove an interesting 
constraint on the origin of the emission-line gas. 

3.3 Density diagnostic 

We note that the [SIII]A 18.71/[SIII]A33.48 line ratio is a density 
diagnostic. In the NGC 1275 East and NGC 4696 pointings we 
have detections of A 18.71 and upper limits for A33.48. Taking 
the ratio of the average of the line fluxes (or upper limits), hav- 
ing used the Ho!-l-[NII] fluxes in the Spitzer apertures to scale the 
fluxes to account for the different spatial sampling as in section lTTI 
gives [SIII]A18.71/ [SIII]A33.48 > 0.6 for NGC 1275 East. Cloudy 
jFerland et ail ll998l) models show that this ratio corresponds to a 
density of > 100 cm^-* . No constraint is placed on the density for 
the NGC 4696 pointing as we are unable to scale the fluxes for the 
aperture effects due to saturation in the Ho!-l-[NII] image. 

3.4 PAH features 

11.3 pLm PAH features are detected in the spectra of NGC 1275 
East and NGC 4696. The flux is strong and comparable to that of 
[Nell] . No detection is made for Perseus positions 2 and 1 1 which 
are at greater distances from the centre of the galaxy than Perseus 
East. The limit on position 2 is less than about one third that of 
[Nell] indicating that either the abundanc e of PAHs o r the e xci- 
tation mechanism is much reduced there. IPeeters et all ( |2004|) ar- 
gue that PAH features are tracers of star formation, pa rticularly of 
B sta rs. Parts of NGC 1275 have an A-type spectrum jMinkowskH 
which indicates some, possibly sporadic, star formation; the 
situation in NGC 4696 is uncertain due to the obvious dust lanes. 



Table 10. The intensity ratio of H2 2 - 1 S(l)/1 - S(l) lines in the three 
regions of NGC 1275. The 2 — 1 S(l) line fluxes are 3f7 upper limits as the 
line is not detected in any region. 



Region 


2-lS(l)/l-0S(l) 


NGC 1275 East 


<0.097 


Position 2 


<0.26 


Position 11 


<0.44 



The correlation between the presence of a PAH feature and a stellar 
continuum in the IRS spectrum in NGC 1275 supports an excita- 
tion mechanism involving the far UV light from moderately young 
stars. 

In a Spitzer study using IRS and Multiband Imaging Photome- 
ter (MIPS) obse rvations of 7 nearby dusty elliptical galaxies includ- 
ing NGC 4696, iKaneda et all ( 120071) find PAH emission from the 
centre and evidence for by dust emission which is more extensive 
than the starlight of the galaxy in the MIPS data. 

3.5 Molecular Hydrogen Lines 

To investigate the molecular hydrogen in more detail we can com- 
bine the mid-infrared pure rotational emission lines with emission 
lines from the I-OS(J) ro-vibrational states visible in the A"— band 
(Section l2.2t 

3.5.1 Excitation process 

Molecular hydrogen can be collisionally excited, or radiatively ex- 
cited by absorbtion of UV photons in the Lyman- Werner bands. 
If the density is low enough (< lO'^cm^-') the molecule will de- 
excite through the ground electronic ro-vibrational states by radia- 
tive de-excitations, whereas if the density is high, the populations 
of the rotational and vibrational states will be redistributed by col- 
lisions before de-excitation to the ground level. The form of the 
observed spectrum depends on the excitation source as well as the 
density and temperature of the gas. At high temperatures and den- 
sities the collision rate of the H2 molecule with other molecules 
and atomic species increases. Collisional de-excitation therefore 
becomes more important than radiative de-excitations as the tem- 
perature and density increase. The total gas density at which the 
collisional de-excitation rate equals the radiative de-excitation rate 
is known as the critic al density and i s a slow function of temper- 
ature (Ster nberg & Da lgarno 1989; M andv & Martinlll993l) . When 
collisional excitation and de-excitation dominates, the gas exists in 
local thermodynamic equilibrium (LTE) and the molecular hydro- 
gen has an ortho-to-para ratio of three. 

3.5.2 Scaling the mid-infrared H2 emission to compare with the 
near-infrared emission 

As discussed in Section|2T|we have assumed that the mid-infrared 
line emission is distributed spatially like the Ha-i-[NII] emission. 
The near-infrared H2 line emission has been shown to c losely fol- 
low t he Pa emission in other brightest cluster galaxies djaffe et all 
|2005^. Additionally, near-infrared H2 emission lines are only de- 
tected in NGC 1 275 in the regions of brightest Ha emission 
dHatch et all2005h . To examine the relationship between the atomic 
hydrogen emission and the near-infrared molecular hydrogen emis- 
sion within NGC 1275 in detail, we have analysed the whole length 
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Figure 9. Pa flux and Pa/1 - S(l) ratio along the UIST long-slit covering 
the NGC 1275 East region. The solid line shows the mean Pa/1 - S(l) of 
1.45. Error bars are at Ifj level. The ratio of atomic-to-molecular hydrogen 
line emission is fairly constant across a large range of Pa emission. 



of the NGC 1275 East HK-hmiA long-slit data. The long-slit data 
were segregated into seven regions in which Pa was dectected. 
Fig.|9]plots the Pa flux and Pa/1 — 0S(1) line ratio against pro- 
jected distance from the galaxy nucleus for all seven regions. This 
figure shows that the Pa/1 — S(l) ratio is almost constant at 1.45 
across a large part of the emission-line region. We therefore as- 
sume that the Pa emission and near-infrared H2 line emission are 
closely related spatially in NGC 1275. We can scale and combine 
the near-infrared and mid-infrared H2 emission lines using the flux 
of the hydrogen recombination lines in the Spitzer and near infra- 
red apertures. 

To obtain Hydrogen line fluxes for the different spectra in dif- 
ferent regions we proceed as follows: For the Spitzer spectra we 
used the Ha-l-[NII] fluxes from the figures given in Table|7]divided 
by two to account for the [Nil] lines ( iHatch et al 2006). For the near 
infra-red spectra we used the Pa flux from Table [9] f or the Perseus 
East region, the Pa flux given by iHatch et al ( l2005h for the Posi- 
tion 2 region and 1 .45 x the flux in the H2 1 — S(l) line (Fig|9]l for 
the Position 1 1 re gion. This corresponds to the 'Horseshoe knot' in 
iHatch et^ ( l2005h which does not have a measured Pa intensity as 
the line was badly affected by the poor atmospheric transmission 
in that part of the spectrum. To scale t he Pa flux to the Ha flux we 
assume Case B recombination theory l lOsterbrock & Ferlandl20d3) 
which implies that Ha/Pa = 8.45. We finally scale the Spitzer Une 
fluxes to match the near infra-red spatial regions with: 



Fs. 



Fnk-(Ha) 
Fs(Ha) 



where: Ess is the resulting scaled Spitzer line flux, Esm is the mea- 
sured Spitzer line flux, Enii-(Ha) is the estimated Ha flux in the near 
infra-red aperture, Es(Ha) is the estimated Ha flux in the Spitzer 
aperture. 



3.5.3 Excitation temperatures, column densities and mass of 
warm molecular hydrogen 

If collisions share the energy between the particles causing the H2 
to be in LTE, the excited states of H2 (N{v,J)) will be populated in 



a thermalised Boltzmann distribution characterised by an excitation 
temperature (Tex) such that: 



8J 



(1) 



where E{v,J) is the upper energy of the {v,J) transition, a{Tex) is 
equal to A'Total/Zr, where A'xotal is the total H2 column density, Zj 
is the H2 partition function, and gy is the statistical weight. N{v,J) 
is the column density of the {v,J) level, which, for optically- thin 
emission can be observationally determined from / (the surface 
brightness [ergs^'cm^^sr^']) through: 



N(vJ) = 



Aui/ic 



(2) 



where A is the rest wa velength of the line and Ayi is the Einstein 
co-efficient taken from iTumer et all l ll977l) . If the H2 level popu- 
lations are completely dominated by coUisional excitation and de- 
excitation then the excitation temperature equals the kinetic tem- 
perature of the gas. 



Population diagrams of In 



gJ 



verses E{v,J)/k\,Tg:-x^ allow 



investigation of whether the levels are thermalised. If the levels are 
thermalised the level populations will lie in a straight line with 
a slope inversely proportional to the excitation temperature. If a 
range of gas temperatures occur along the line-of-sight, the points 
will lie on a smooth curve with the lower energy levels lying on 
a steeper slope (lower temperature) than the higher energy levels 
(higher temperature). This is because the higher energy levels are 
preferentially populated at higher temperatures. 

Population diagrams for the NGC 1275 East, position 2 and 
position 1 1 regions are presented in Eigs.[T0l[TT]and[T2] As the pop- 
ulation of the states lie on a smooth-curve, these diagrams clearly 
show that the molecular hydrogen exists in LTE at a range of tem- 
peratures. Supporting evidence that the H2 is thermally excited 
comes from the ratio of the H2 2— 1S(1)/1— 0S(1) line intensity 
ratio. 3cT upper limits for the 2—1 S(l) emission lines were mea- 
sured from the position 2, 1 1 and the NGC 1275 East regions and 
the limits on the 2— 1S(1)/1— 0S(1) ratios are given in Table 1 101 
The measured ratio is si gnificantly le ss than the pure radiative de- 
excitation value of 0.53 (Mouri 1994), so the gas is predominately 
de-excited through collisions, perhaps with a small non-thermal 
component. This implies that the gas must have a density of greater 
than 10"^ cm"^. 

To gain a simple picture of the structure of the molecular gas 
we have measured the rotational temperature from the — 0S(1) 
and — S(2) lines, and a ro- vibrational temperature from the 1-0 
lines. In the NGC 1275 East HK-hand spectrum we measured the 
I — 0S(7) line, and therefore we measure two ro-vibrational ex- 
citation temperatures in this region: one from levels with an up- 
per energy of < 8000 K and one from levels with an upper energy 
of > 8000 K. These excitation temperatures are given in Table [TT] 
The pure-rotational temperatures of the gas are 330 — 370 K for 
all regions. The ro-vibrational lines reveal a layer of hotter gas at 
1700 - 2000 K and in NGC 1275 East there is an additional layer of 
even hotter gas at ~ 2600 K. Such high population temperatures are 
often int erpreted as due to no n-thermal excit ation, for instance b y 
starlight jvan Dishoeckl2004h or cosmic rays toalgamo et all 199^ . 
In the remainder of this paper we will refer to the high v,J popula- 
tions as an indication of warm H2, without meaning to specify what 
actually causes these level populations. 

The total H2 column density of the molecular hydrogen can 
be calculated through rearranging equation[T]and assuming that all 
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Table 11. Excitation tempertures, column densities and masses of molecular Hydrogen. The observed mass is calculated directly from the observed flux in the 
relevant aperture. The scaled mass is the mass scaled to the region in which the ro-vibrational lines were observed. Error bars are at the la level. 
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Figure 10. Population diagram of the NGC 1275 East region with Iff eiTor 
bars. The datapoints at £" = 1015 K and 1682 K are measured by Spitzer and 
the error bars include the uncertainty in scaling the mid-infrared emission 
to the near-infrared emission using hydrogen recombination lines. 
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Figure 11. Population diagram of NGC 1275 position 2 with Iff error bars. 
The datapoints at£' = 1015K and 1682 K are measured by Spitzer and the 
error bars include the uncertainty in scaling the mid-infrared emission to 
the near-infrared emission using hydrogen recombination lines. 



the molecular hydrogen exists in LTE at the excitation temperature 
derived from the population diagrams. 

N{v,J)Z{T) 



gje 



(3) 



where Z{T) is the partition function. 

As the molecular hydrogen exists at multiple temperatures 
we have assumed the gas consists of two (or three in the case of 
NGC 1275 East) separate populations, each consisting of a Boltz- 
mann distribution at the excitation temperatures derived from the 
population diagrams. A least squares fit was performed with the to- 
tal column densities of each Boltzmann distribution as the two (or 
three) free parameters. Column densities at each excitation temper- 
ature for the three regions are given in Table [TT] 



The mass of warm molecular hydrogen within the Spitzer 
aperture can be derived from a single line luminosity and the ex- 
citation temperature. The total molecular hydrogen mass Mxotal is 
derived from: 



MTotal =MH,ntotal, 



(4) 



where Me, is the mass of a hydrogen molecule, and njotal is the total 
number of H2 molecules, njotal can be derived from the luminosity 
from the (v,J) state through: 



Utotal ■ 



(5) 



where L{v,J) is the luminosity in the {v,J) line, Zj^^ is the partition 
function, A^i is the Einstein co-efficient, and V is the frequency of 
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Figure 12. Population diagram of NGC 1275 position 11 with Iff eiTor 
bars. The datapoints at £" = 1015 K and 1682 K are measured by Spitzer, 
but the error bars do not include any uncertainty associated with scaling the 
mid-infrared emission lines. 



the line. The masses of the ---^350 K material in each of the Spitzer 
apertures were calculated from the 0-0 S(l) line. The masses of the 
warmer gas (1700-2000 K) were measured using the l-OS(l) line. 
A check was made using the other detected lines 1-OS (2), 1-OS (3), 
1-OS (4) and the masses derived from these lines were within 10 per 
cent of the massess derived from the l-OS(l) line. Finally the mass 
of the 2580K material seen in the NGC 1275 East region was mea- 
sured from the 1-0 S(7) line. The masses of molecular hydrogen 
derived from the different temperature lines are given in Table [TT] 
as the "Observed mass of H2". In the rows labelled "Scaled mass of 
H2" we list the masses of molecular hydrogen at each temperature 
scaled to the region in which the ro-virational lines were measured; 
the Spitzer masses were scaled down by the ratio of Ha fluxes in 
the relative apertures, as described in Section [3.5.2l The majority 
of the molecular gas exits at lower temperatures with less than 0.5 
per cent of the total mass being at the higher temperatures. 

Without the measurement of the 0-0 S(0) line, it is not possi- 
ble to derive the temperature and mass of of the coldest component 
of the warm H2 gas. Therefore we use the 3cT upper limits on the 
— S(0) flux to obtain lower limits on the temperature and upper 
limits on the mass of the coldest component of the warm molecular 
hydrogen. These limits are displayed for the NGC 1275 pointings in 
Tablell2l We cannot carry out the scaling for the NGC 4696 point- 
ing as the region of the Ha-l-[NII] map covered by the LH aperture 
is affected by saturation. The minimum allowed temperature of the 
coolest component is 150-180 K, and approximately 10 times more 
mass may exist at this cool temperature than in the warmer molec- 
ular hydrogen. 



3.5.4 Correlation with Optical Emission Lines 

In Fig. [T3]we plot the molecular hydrogen 0-0 S(l) flux against 
the Ha-l-[NII] flux for all the NGC 1275 and NGC 4696 pointings. 
The flux in the Spitzer 0-0 S(l) molecular hydrogen line seems to 
correlate well with the flux in the optical Ha-l-[NII] lines. The fit is 
a straight line constrained to go though the origin and shows that the 
flux in the O-OS(l) line is 0.03 x the flux in the Ha-l-[NII] complex. 




NII)xlO-'-'ergcm-V' 



Figure 13. Flux in 0-0 S(l) Une plotted against the flux in the Ha-H[NII] 
complex for all the NGC 1275 and NGC 4695 pointings. The solid line is 
a straight-line fit in linear coordinate space, constrained to pass though the 
origin. It has a slope of 0.03. 



3.5.5 Pressure balance between the molecular, atomic and X-ray 
emitting gas 

There is a serious problem implied by the apparent thermal dis- 
tribution of the levels in the molecular hydrogen. This is that 
the currently used coUisional rates for the r o-vibrational states re- 
quire that the density exceeds 10"*^^ cm^^ l lStemberg & Dalgamd 
1989; Mandv & Martin 199^). If the temperature is few x 1000 K 
then the pressure is 10^^*^ cm^"* K. This exceeds the ambient pres- 
sure, as determined from the X-ray measurements of the hot gas 
jSanders et aill2005h . by an order of magnitude or more. This does 
not appear to be a stable situation for such thin, apparently long- 
lived, filaments. 

Recent studies have brought the H2 collision database into 
question. The Orion Bar is t he closest and bes t studied H+ / H" 
/ H2 interface ( IO'Deli200lh . lAllers et all ( |2005|) found that H2 ro- 
tational lines implied level populations that were close to a thermal 
distribution. They noted that the density within the Bar is too low 
to thermalize the levels involved if the current H2 collision rates are 
correct. They suggest that current H - H2 vibrational de-excitation 
rates should be increased by nearly two orders of magnitude (their 
Table 6). If this were the case then H2 level populations could reach 
a thermal distribution at the low densities implied by the pressure 
of the hot gas. Our data are further evidence that the de-excitation 
rates may require a large correction. 



3.6 Heating Mechanisms 

A detailed discussion of the heating mechanism of the emission- 
line gas in central galaxies in cooling core clusters is beyond the 
scope of this paper. In the nuclear regions the situation can be 
complex with several different mechanisms at work. These in- 
clude photo-ionizat ion from an active nucelus (for NGC 1275) (eg 
Wilmaneti3l2005l ). photoi onization b y stars (eg Ijohnstone et all 



19871; 



Allerj|l99 5: Wilma n et all l2002h and possibly shocks (eg 



iJaffe et alll200ll ; IWilman et alll2002l)" 



We expect that the off-nuclear regions which we have targeted 
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Table 12. Limits on the temperature and molecular hydrogen masses from the O-OS(O) lines. The observed mass limits are calculated directly from the upper 
limit to the observed flux in the relevant aperture. The scaled mass is the mass scaled to the region in which the ro-vibrational hues were observed. 



Region Temperature Observed mass of H2 Scaled mass of H2 

NGC 1275 East >180K <2.9xlO^M0 < 3.2 x 10^ Mq 

NGC 1275 Position 2 > 155 K <1.Ox1O'M0 <3.8x10'Mq 

NGC 1275 Position 11 > 150 K < 5.6 x lO*" Mq 0.2x10' M,.;, 



in the observations presented in this paper will be more represen- 
tative of the general emission-line nebula and also simpler to un- 
derstand. We note that the apparently thermalized nature of the 
molecular hydrogen lines argues against their formation in a clas- 
sical photodissociation region and that the relation between the 
Ha+[NII] flux and the flux in the molecular hydrogen lines pre- 
sented in Fig. [T3] could suggest that shocks are important. How- 
ever, it is not clear that there would be enough cloud-cloud colli- 
sions in these regions to drive this mechanism. Shocks also need 
to be constrained to a fairly narrow range of velocity, fast enough 
to produc e the optic al emission lines, but not so fast that dust is 
destroyed. IJaffe et al ( .2001.) have previously discussed the require- 
ment for the fine-tuning of shock velocities in trying to explain the 
ro-vibrational molecular lines in a sample of cooling core central 
cluster galaxies. 

In a future paper we will explore these difficulties in more de- 
tail and examine the possibility of excitation via cosmic rays which 
leak out of the ra dio lobes and ghost bubbles (eg lBirzan et all2004l : 
iDunn et afcOOSi) . 



energy flow in the filaments. However the total molecular hydro- 
gen emission, including the near IR ro-vibrational lines is compa- 
rable to the emission from Ha. (Ha is about one tenth of the total 
line emission: the major emitter is expected to be Lya.) Neverthe- 
less they reveal more about the low temperature core to the fila- 
ments which is dominated in mass by H2. The total mass in H2 for 
NGC 1275, obtained from the CO emission using a standard Galac- 
tic conversion ratio, is '-^ 4 x lO'*' Mm (Salome et al 2006), most of 
it near the centre of the galaxy. The mass in H2 at 300 K is about 
6x 10^ M,i„ assuming that the masses found for the outer filaments 
in this paper scale linearly with Ha emission, and using the total 
Ha flux of Conselice et al (2001). The total mass of Ha emitting 
gas is smaller at about 3x10^ Mq and the mass of surrounding soft 
X-ray emitting gas, with temperature ~ 7 x 10^ K, is about 10^ Mq 
(Fabian et al 2006). 

Further observations are required to test whether all filaments 
have similar composition and properties, but they do appear to be 
predominantly molecular. 
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4 CONCLUSION 

Using high resolution mid-infrared spectra from the Spitzer 
Space Telescope we have detected off-nuclear emission lines in 
NGC 1275 and NGC 4696. The lines arise from pure rotational 
transitions of the Hydrogen molecule, atomic fine-structure tran- 
sitions from Ne+, Ne++, S++ and Si+ and the 11.3 micron PAH 
feature. 

The relative intensities of the molecular hydrogen lines are 
consistent with collisional excitation and probe a new region of 
temperature space in the filaments at 300-400 K. Pressure equilib- 
rium with the surrounding intracluster medium requires a substan - 
tial revision to H2 collision rates, as noted by I AUers et all ( l2005h . 
The outer filaments around NGC 1275 therefore have molecular 
gas at 50 K (CO, from P. Salome et al, i n preparation), 30 0-400 K 
(H2, this paper) and 2000-3000 K (H^.lHatch et alll2005l) , atomic 
as at seve ral 1000 K (Ha, e.g. lConselice et alll20011) "and ~ 10^ K 
Fabian et a l 2003b), embedded in the hot intracluster gas at about 

5 X 10^ K. Detection of OVI from more central regions with FUSE 
teregman et all2003) suggests that gas at ~ 3 x 10^ K is also likely 
to be present. 

The discovery of the [Nelll] lines is unexpected based on the 
weakness of the optical [OIII] lines and suggests that the elec- 
tron temperature in the filamentary medium where these lines are 
formed may be very low, possibly due to enhanced metallicity. The 
metal abundance of the central intracluster medium in cooling clus- 
ter cores is enhanced by Type la, and possibly Type II supemovae 
(see Sanders & Fabian 2006 and references therein), so the fila- 
ments are enriched if they orginate from cooling of that gas. 

Collectively the mid-IR emission lines radiate about one quar- 
ter of the flux in Ha, so they play a relatively minor role in the 
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